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ENGINE CONTROL USING TORQUE ESTIMATION 

Inventors: Giorgio Rizzoni 
Byungho Lee 
Yann Guezennec 
Ahmed Soliman 
Mauro Cavaelletti 
James Waters 

This application claims the benefit of U.S. Provisional Patent Application No. 
60/273,423 entitled ENGINE CONTROL USING TORQUE ESTIMATION and filed 
March 5, 2001. 

TECHNICAL FIELD 

The present invention relates to systems and methods for engine control. In 
particular, the present invention relates to a system and method for engine control using 
stochastic and frequency analysis torque estimation techniques. 

BACKGROUND AND SUMMARY OF THE INVENTION 

In recent years, the increasing interest and requirements for improved engine 
diagnostics and control has led to the implementation of several different sensing and 
signal processing technologies. In order to optimize the performance and emission of 
an engine, detailed and specified knowledge of the combustion process inside the 
engine cylinder is required. In that sense, the torque generated by each combustion 
event in an IC engine is one of the most important variables related to the combustion 
process and engine performance. 



In-cylinder pressure and engine torque have been recognized as fundamental 
performance variables in internal combustion engines for many years now. 
Conventionally, the in-cylinder pressure has been directly measured using in-cylinder 
pressure transducers in a laboratory environment. Then, the indicated torque has been 
calculated from the measured in-cylinder pressure based on the engine geometry while 
the net engine torque has been obtained considering the torque losses. However, such 
direct measurements using conventional pressure sensors inside engine combustion 
chambers are not only very expensive but also not reliable for production engines. For 
this reason, practical applications based on these fundamental performance variables in 
commercially produced vehicles have not been established yet. Therefore, instead of 
employing the expensive yet not reliable conventional approach, there is a need for 
different approaches of obtaining and using such performance variables by estimating 
the net cylinder torque resulting from each combustion event while utilizing pre-existing 
sensors and easily accessible engine state variables, such as the instantaneous 
angular position and velocity of the crankshaft. This approach enhances the on-board 
and real-time estimations of engine state variables such as instantaneous torque in 
each individual cylinder and bring out many possible event-based applications for 
electronic throttle control, cylinder deactivation control, transmission shift control, misfire 
detection, and general-purpose condition monitoring and diagnostics [1-3]. 

The crankshaft of an IC engine is subjected to complex forces and torque 
excitations created by the combustion process from each cylinder. These torque 
excitations cause the engine crankshaft to rotate at a certain angular velocity. The 
resulting angular speed of engine crankshaft consists of a slowly varying mean 



component and a quickly varying fluctuating component around the mean value, caused 
by the combustion events in each individual cylinder [4]. Outcome of the torque 
estimation approaches strongly relies on the ability to correlate the characteristics of the 
crankshaft angular position, speed, and its fluctuations to the characteristics of actual 
cylinder torque [3] and [4]. Over the past years, this torque estimation problem has 
been investigated by numerous researchers explicitly or implicitly, inverting an engine 
dynamic model of various complexities. Those researchers have successfully 
developed and validated the dynamic models describing the cylinder torque to the 
crankshaft angular velocity dynamics in internal combustion engines. 

One of the earliest strategies targeted at developing the engine and crankshaft 
dynamic model allowing the speed-based torque estimation was carried out by Rizzoni, 
who introduced the possibility of accurately estimating the mean indicated torque by a 
two-step procedure [4]. It consists of first deconvolving the measured crankshaft 
angular velocity through the rotational dynamics of the engine to obtain the net engine 
torque which accelerates the crankshaft, and then of converting this net torque to 
indicated torque through a correction for the inertia torque component, caused by the 
reciprocating motion of crank-slider mechanism, and for piston/ring friction losses. 
Another strategy was introduced focusing on reconstructing the instantaneous as well 
as average engine torque based on the frequency-domain deconvolution method [3]. 
However, this method required pre-computation of the frequency response functions 
relating crankshaft speed to indicated torque in the frequency-domain and storing their 
inverses in a mapping format, which has difficulties of determining the frequency 
functions experimentally. An approach bypassing this difficulty was proposed by 



Srinivasan et al. using the repetitive estimators [5]. Further studies of the speed-based 
torque estimation was continued by Kao and Moskwa, and Rizzoni et al. through the 
use of nonlinear observers, particularly sliding mode observers [6] and [7]. This method 
of the nonlinear observer was desirable for variable speed applications since a wide 
range of operating conditions required the non-linearity of the models. Other torque 
estimation efforts involving an observer were based on the use of the unknown input 
observer by Rizzoni et al. [8-10]. This method was, however, only applicable to 
constant speed (or near constant speed) engines. One of the most recent research 
efforts aimed at the individual cylinder pressure and torque estimations was based on 
the stochastic approach by Guezennec and Gyan [1] and [11]. This approach permitted 
estimations of the instantaneous in-cylinder pressure accurately without any significant 
computational requirement based on the correlations between in-cylinder pressure and 
crankshaft speed variations. 

Even though all these approaches described previously were successful over the 
past years, most of them were not feasible for the on-board real-time estimation and 
control in mass-production engines. In other words, these approaches can only be 
practically implemented in a post-processing phase because they must involve either 
highly resolved measurements of the crankshaft speed or significant amounts of 
computational requirements. The present invention, however, presents a practical and 
applicable way of implementing the speed-based torque estimation technique on a 
production engine in order to develop a methodology and algorithm extracting the in- 
cylinder pressure and indicated torque information from a less resolved/sampled 
crankshaft speed measurement for the purpose of real-time estimation and engine 



control in production vehicles. Two different approaches have been implemented, 
namely "Stochastic Estimation Technique" and "Frequency-Domain Analysis," to 
estimate the instantaneous indicated torque (as well as in-cylinder pressure) in real time 
based on the crankshaft speed fluctuation measurement. An overview of both 
techniques is presented. Next, their implementations on an in-line four-cylinder spark- 
ignition engine are presented under a wide range of engine operating conditions such 
as engine speed and load. Then, validations of the robustness of these techniques are 
presented through the real-time estimation of indicated torque during the actual engine 
9 operations, demonstrating that these methods have very high potential for event-based 

:|! engine controls and diagnostics in mass-production engines. 
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f BRIEF DESCRIPTION OF THE DRAWINGS 

yj Figure 1 is a Simplified SISO Model for Engine Dynamics for an example 

W embodiment of the present invention; 

m 

f# Figure 2 shows Basis Variables for Pressure Estimation for an example 

embodiment of the present invention; 

Figure 3 shows an In-Cylinder Pressure Estimation at Speed of 2000 RPM and 
Load Torque of 30 Ibf-ft for an example embodiment of the present invention; 

Figure 4 shows an In-Cylinder Pressure Estimation for an example embodiment 
of the present invention; 

Figure 5 shows Indicated Torque Estimation for Each Cylinder for an example 
embodiment of the present invention; 
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Figure 6 shows Indicated Torque Estimation for All Cylinders for an example 
embodiment of the present invention; 

Figure 7 shows Indicated Torque Estimation for Each Cylinder for an example 
embodiment of the present invention; 

Figure 8 shows Indicated Torque Estimation for All Cylinders for an example 
embodiment of the present invention; 

Figure 9 shows Cycle-Averaged Indicated Torque Estimation for an example 
embodiment of the present invention; 

Figure 10 shows Average R.M.S. Errors for Various Cases for an example 
embodiment of the present invention; 

Figure 11 shows Spatial Spectra for Indicated Torque for an example 
embodiment of the present invention; 

Figure 12 shows Spatial Spectra for Speed Fluctuation for an example 
embodiment of the present invention; 

Figure 13 shows Coherence Function for Crankshaft Speed Fluctuations and 
Indicated Torque for an example embodiment of the present invention; 

Figure 14 shows Average Indicated Torque vs. Approximated R.M.S. of Torque 
Fluctuations for an example embodiment of the present invention; 

Figure 15 shows Indicated Torque Estimation at 2000 RPM and 53 N-m Load 
Torque for an example embodiment of the present invention; 

Figure 16 shows Coefficient Estimation at All Operating Points for an example 
embodiment of the present invention; 



Figure 17 shows Indicated Torque Estimation of Each Cylinder for an example 
embodiment of the present invention; 

Figure 18 shows Indicated Torque Estimation of All Cylinders for an example 
embodiment of the present invention; 

Figure 19 shows R.M.S. Error for Various Cases for an example embodiment of 
the present invention; 

Figure 20 shows Real-Time Estimation of Individual Cylinder Torque for an 

M> example embodiment of the present invention; 

CI 

Figure 21 shows Actual Value of Indicated Torque from Acquired Data for an 

$ 

y example embodiment of the present invention; 

]JJ Figure 22 shows Real-Time Estimation of Summation of Indicated Torque for an 

,n 

ip example embodiment of the present invention; and 

Mi 

ffl Figure 23 show Actual Value for Sum of Indicated Torque from Acquired Data. 
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DETAILED DESCRIPTION OF EXAMPLE EMBODIMENTS 

STOCHASTIC ESTIMATION TECHNIQUE 
This technique is based on a signal processing method, herein referred to as the 
"Stochastic Estimation Method," which allows extraction of reliable estimates based on 
the method of least square fittings from a set of variables which are statistically 
correlated (linearly or otherwise). The procedure originates from the signal processing 
field, and it has been used in a variety of contexts over the past years, particularly in the 
field of turbulence [1]. It has been primarily used for estimating conditional averages 
from unconditional statistics, namely, cross-correlation functions. The main advantage 
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of this methodology compared to others is that all complexities of the actual physical 
system are self-extracted from the data in the form of first, second, or higher correlation 
functions. Once the correlation models are determined, the estimation procedure 
reduces to a simple evaluation of polynomial forms based on the measurements. 
Consequently, the estimation can be achieved in real time with very few computational 
operations. The stochastic estimation methodology may be used in order to achieve the 
estimation of in-cylinder pressure and indicated torque based on the crankshaft speed 
measurements. 

ft 

0 A given set of variables of x 1} x 2 , x 3 , and x 4 may be statistically correlated with 

m . 

111 another variable of y. Each variable has N number of realizations or measurements. 

0 

* Then, a polynomial equation to express y in terms of x 1 though x 4 can be written as 

\wh ■ 

U y estimate = a 0 + <*l X l + a 2 X 2 + a 3 X S + a 4 X 4 ( 1 ) 

yl 

& where ao to a 4 are the polynomial coefficients. Applying the least mean squares gives 

m .. 

SI the expression of an error between the true value of y (y tru e) and estimated value of y 

(yestimate) SUCh as 

N 

8 = J] OW* ~ y estimate* ) (^) 

where ^ is the estimation error, and N is the total number of realizations. Then, the 
polynomial coefficients in Eq. (1), a 0 through a 4 , must be determined so that Eq. (1) 
estimates the variable y as best as possible based on the statistical sample of N 
realizations. This best estimation corresponds to minimizing the error term s over all 
realizations, which leads to taking the partial derivatives of the error in Eq. (2) with 
respect to each of the coefficients and then setting them equal to zero. This procedure 
results in the following set of equations. 
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a 0 Yl + aj-Lx IJt +a 2 1.x 2jc + a 3 T.x 3Jl + a 4 l,x 4jl 

+a 7 l4 +o^uPhx =5>w>W 
flbl^ +^2^ =2>at)W 

+<%Z%% +^2x^t% +<^Z^* = £fyt>W 
Taking an average over all realizations for each equation then converting them into a 

matrix form gives the following final format. 

" 0 {*,) fe> (*<) ' 

( X j) (xj) (x 1 x 2 ){x 1 x 3 ){x 1 x 4 ) 

(x 2 ) (x t x 2 ) (x 2 2 ) {x 2 x 3 )(x 2 x 4 ) 

(x 3 ) (x,x 3 )(x 2 x 3 ) (x 2 3 ) (x 3 x 4 ) 

_ (*<) (xjX 4 ){x 2 x 4 }(x 3 x 4 ) (x 2 ) 

where <> denotes averaging over all realizations. After the cross-correlation matrices 
have been constructed based on all the available N realizations as shown in Eq. (3) 
above, the set of polynomial coefficients, a 0 through a 4 , can be determined once for all. 
Then, the variable yean be estimated using Eq. (1) during the estimation phase without 
any significant computational requirement. For the implementations of this technique on 
IC engines, it is necessary to obtain quantitative representations of the in-cylinder 
combustion events, such as in-cylinder pressure and indicated torque, based on the 
given measurements of the crankshaft rotational dynamics (position, speed, and 
acceleration). Therefore, cross-correlation functions may be built as shown in Eqs. (1) 
and (3) between the quantities to estimate (in-cylinder pressure or indicated torque) and 
the quantities measured (or combinations of those quantities). 

FREQUENCY ANALYSIS TECHNIQUE 
One of the main advantages of using the frequency domain technique is that the 
accuracy of the estimation can be improved by performing the operation in the 



"<«.> 

fa) 
(*) 
(*,) 



{ytrue) 

{y*m*i) 
(y^Xs) 



(3) 
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frequency domain rather than in the time or crank angle domain, considering only a few 
frequency components of the measured crankshaft speed signals [3]. This 
reconstruction technique is feasible mainly due to the intrinsically periodic nature of the 
engine process, which leads to the use of Fourier Transform as a tool of performing the 
crankshaft speed deconvolution through the engine crankshaft dynamics. The 
computation in the frequency domain, employing the Discrete Fourier Transform, 
effectively acts as a comb filter on the speed signal and preserves the desired 
information, which is strictly synchronous with the engine firing frequency [3]. This 
frequency domain deconvolution is very effective mainly because it reduces the process 
to an algebraic operation and the dynamic model representing the rotating assembly 
needs to be known only at the frequencies that are harmonically related to the firing 
frequency [4]. 

In order to perform the speed-based torque estimation using the frequency 
approach, the engine crankshaft dynamics are considered as a SISO (Single-Input & 
Single-Output) model, as described in Fig. (1). 

Within Fig. (1), the indicated torque (denoted by T(Q)) is considered as an input 
to the engine dynamic system (denoted by H(9)), and the crankshaft speed (denoted by 
Q(9)) is considered as a system output resulting from the torque generated by the 
engine. Because those signals are acquired in the crank angle domain as denoted, the 
Fourier Transform generates the spatial spectrum. The relationship between the 
indicated torque and crankshaft speed in the spatial frequency domain can be described 
as shown in Eq. (4) below 

T l (M = a(jA k )H-'(jA k ) (4) 
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where j is the imaginary part, fa is the angular frequency (/c^ order of rotation), rfiAk) 
and £2(jXk) are the Fourier Transforms for the indicated torque and crankshaft speed 
respectively, evaluated at a frequency of A k , and HQty is the engine frequency 
response function evaluated at that frequency. Therefore, the frequency response 
function H is obtained at each of the first few harmonics of the engine firing frequency 
through either experimental data or theoretical models. Then, computing the Discrete 
Fourier Transform of the crankshaft speed {£2Q2,$) at each of the selected harmonics 
allows us to evaluate the indicated torque in the frequency domain (n(j^) at each 
harmonic using Eq. (4). Finally, t,QA.^ can be converted into the crank angle domain 
using the Inverse Discrete Fourier Transform at each of the harmonics in order to obtain 
the estimation of the indicated torque. To implement this approach on IC engines in 
real-time, the first few harmonics of the firing frequency within the signals contain 
enough information to represent the actual engine behavior between the crankshaft 
speed and indicated torque of the simplified SISO engine dynamics model described in 
Fig.(1)[4]. 

EXPERIMENTAL DATA 
In order to validate and implement the approaches described previously, the 
estimation techniques were applied to a set of experimental data acquired from a 2.4L, 
DOHC, in-line four, spark-ignited, passenger car engine manufactured by General 
Motors. The main characteristics of the engine are described in Table (1) below. 
Results from this data set are provided. The experimental data sets consist of various 
measurements, listed in Table (2), with an angular resolution of 1° of crank angle (720 
data points per engine cycle) and 100 consecutive engine cycles for each 
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measurement. Each data set was acquired under a wide range of engine operating 
conditions for various engine speed and load, as shown in Table (3). 
Table 1 . Characteristics of Engine 



Engine Type 


1-4, spark 
ignited, DOHC 


Bore 


90 mm 


Stroke 


94 mm 


Connecting Rod Length 


145.5 mm 


Displacement Volume 


2.4 liter 


Number of Valve 


4 per cylinder 


Compression Ratio 


9.7 



Table 2. List of Measured Data 



TDC of Cylinder #1 


Intake Air Flow Rate 


Each Cylinder 
Pressure 


Load Torque 


Crankshaft Speed 


Intake Air 
Temperature 


Intake Manifold 
Pressure 


Exhaust Gas 
Temperature 


Air/Fuel Ratio 


Engine Oil 
Temperature 


Spark Ignition Timing 


Coolant Temperature 


Fuel Injection Timing 


Throttle Position 



Table 3. Various Engine Operating Conditions 



Load Torque 
[Ibt-ft] 


Engine Speed [RPM] 

(With an Increment of 500 
RPM) 


10 


1000 to 5000 RPM 


30 


1000 to 5000 RPM 


50 


1500 to 5000 RPM 


70 


2000 to 5000 RPM 
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90 



2000 to 5000 RPM 



TORQUE ESTIMATION USING STOCHASTIC ANALYSIS METHOD 
A direct application of this methodology on the speed-based torque estimation is 
described. There are two separate approaches to estimate the indicated torque based 
on the crankshaft speed fluctuations. The first approach consists of estimating the in- 
cylinder combustion pressure then calculating the indicated torque based on the 
estimated pressure and the engine geometry. The other approach consists of directly 
U estimating the indicated torque from the crankshaft speed fluctuation measurement. 

In any case of estimation approaches, the estimation model function (referred as 
II the basis function) consists mainly of three primary variables representing the 



m 

m 



W crankshaft dynamics such as crankshaft position, speed, and acceleration. A function 
related to the crankshaft angular position is included instead of crank angle itself in the 
basis function because the angular position is clearly cyclic with a period of 4n thus 
introduces a discontinuity at every engine cycle. Because the mathematical foundations 
of the stochastic technique are continuous in nature, this discontinuity leads to 
undesirable mathematical errors. Consequently, a function that is mathematically 
related to the crankshaft position but more closely related to the behaviors of in-cylinder 
pressure or indicated torque is more appropriate. Because the compression and 
expansion strokes, excluding the combustion event, can be considered as polytropic, 
the in-cylinder pressure roughly follows pV* = constant [12]. Because the volume of a 
cylinder for a given engine can be easily obtained from the given engine geometry and 
measured crank angle, a position function U can be considered to be directly 
proportional to V k during the compression and expansion strokes, and constant 
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elsewhere in order to represent the position of the crankshaft [1] and [11]. Such 
function has a high level of correlation with the measured in-cylinder pressure or with 
the measured indicated torque since it effectively represents the motored pressure or 
motored torque information. For the crankshaft speed signal, the relevant signal is the 
crankshaft velocity signal fluctuating around its mean value. Therefore, the general 
correlation function for estimating the in-cylinder pressure or indicated torque can be 
written as a function of the position function f ft angular speed fluctuation 6, and angular 
acceleration 6, as shown below. 

Estimated Value = f[/ 0 $d) (5) 
ESTIMATION OF IN-CYLINDER PRESSURE 
After the in-cylinder combustion pressure is estimated based on the crankshaft 
speed measurement, the indicated torque is then calculated accordingly based on the 
estimated in-cylinder pressure and the given engine geometry. The estimation model 
function (basis function) may be set to be the following first-order non-linear model as 
shown in Eq. (6) in order to first estimate the in-cylinder pressure. 

^estimate 

= a„ +a t f 0 + a 2 f e (Ka 3 f e $¥a 4 (k (6) 

The stochastic estimation approach requires building the cross-correlation 
functions between the estimation quantity (in-cylinder pressure) and the measured 
quantities (three basic variables as well as their cross-terms as shown in Eq. (6)). The 
coefficients, a 0 through a 4 , can be obtained by minimizing the mean square difference 
between the measured pressure and the estimated pressure as shown in Eq. (7). 
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As described earlier in Eqs. (2) and (3), taking the partial derivatives with respect 
to each of the coefficients and setting the result equal to zero gives the following cross- 
correlation matrix system to solve. 

<;> <f e > <f e $> <f $ $> <M> 

<fe> <fo> <//* 

<f&> </;$> <rfi > 
<fA<fi<& <fi&><fjfa> 
<M><f g M>< > < > 

In Eq. (8), the various terms in the matrix represent the cross-correlations among 
the measured basis variables while the right side of the equation represents the cross- 
correlations between the measured in-cylinder pressure and the measured basis 
variables. These non-linear cross-correlations are pre-computed based on all available 
data at a certain engine operating condition, then the five coefficients are computed 
once for all (cycles and cylinders) at that operating point. Once the coefficients as well 
as these correlation functions are determined and proper processing has been carried 
out, the estimation procedure reduces down to the simple evaluation of a multivariate 
polynomial form based on the measurements. Therefore, during the estimation phase 
the instantaneous value of the five measured basis variables are used to evaluate the 
simple polynomial equation as shown in Eq. (6) for the desired estimation. Therefore, 
the computational requirements can become very minimal in this approach, and the 
estimation can be achieved in real time with a few computational operations. 

Referring to Fig. (2), Fig. (2) represents each of the prescribed basis variables 
including the in-cylinder combustion pressure position function U Based on these 
variables, the in-cylinder pressure was estimated using the basis function described in 
Eq. (6) and the cross-correlation described in Eq. (8). Referring to Fig. (3), Fig. (3) 
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(8) 
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represents the estimated in-cylinder pressure trace in comparison with the measured 
trace at a certain engine operating point 

Referring to Fig. (3), the in-cylinder pressure estimation closely follows the 
actually measured pressure trace for each of the cylinders with only minor errors. 
Based on the estimated pressure and the given engine geometry shown in Table (1), 
the individual cylinder indicated torque and summation of the individual cylinder torque 
can be calculated as well [12]. 

However, this estimation is based on the resolution of 360 per crankshaft rotation 
(every 1° of crank angle), which would require a substantial computation power for the 
real-time estimation purpose. For this reason, using fewer resolved measurements, 
such as 36 and 60 resolutions, may allow this technique to be feasible for the real-time 
estimation and control application. Figure (3) represents the in-cylinder pressure 
estimation based on the 36 resolutions (every 10° of crank angle). 

Referring to Fig. (4), using fewer sampled measurements during the computation 
can also provide a successful in-cylinder pressure estimation just as using the full 360 
resolutions can. Based on this pressure estimation and the given engine geometry, the 
individual cylinder indicated torque and summation of the individual cylinder torque were 
calculated and are shown in Figs. (5) and (6), respectively. 

In order to compare the estimation accuracy of different resolutions and possibly 
different estimation models in the later analysis, an error function was defined as the 
root mean square (R.M.S.) error between the measured pressure and estimated 
pressure. Then, this R.M.S. error was normalized by the peak pressure averaged over 
all cylinders and cycles, as shown in Eq. (9) below. 
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Normalized R.M.S. Error = 



V v^esM Pmeasj ) 



(9) 



Table (4) illustrates this estimation error for each of the estimations and number 
of resolutions accounted in the computation. Note that the values are averages over all 
engine operating conditions. 
Table 4. Normalized R.M.S. Errors for Various Cases 



Estimation 
Type 


Number of Resolutions 


360 


60 


36 


Indicated 


Pressure 


2.694 % 


5.063 % 


3.494 % 


Indicated 
Torque 


Individual 
Cylinder 


3.394 % 


5.810 % 


4.313 % 


All 
Cylinder 


6.159% 


7.603 % 


6.814 % 



ESTIMATION OF INDICATED TORQUE 
The indicated torque is estimated directly from the crankshaft speed 
measurements, replacing the two steps procedure of first estimating the in-cylinder 
pressure and secondly calculating the indicated torque accordingly.' There are two 
different parts of achieving the indicated torque estimation in this approach. The first 
part is to estimate the individual cylinder torque for each cylinder then calculate their 
summations whereas the other part is to directly estimate the summation of individual 
cylinder torque. 

Basis Function Selection - Various basis functions are investigated in order to 
determine the best form of the estimation model for the indicated torque estimation in 
real-time. 

Table 5. Various Basis Functions 



Function 
Number 


Basis Function 


1 
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2 


T aMaK =a„+a,f l) +a 2 <k-a 3 &a 4 $ 


3 


T M =a 0 + a,f, + a 2 f e A a s fA affi 


4 




5 


T, SWM = a 0 + aj e + a ? $f a^k aj e 2 + a$ + a 6 0 2 


6 


T eimw =a„+ a,f, + a 2 fjh a 3 f e <& a$ + a s M- a 6 $ 




Tesmx: = a 0+ "iff + «2<^f a 3^ «j} 


7 


+a s f e <&-a 6 f,Aa r $ +a s $$-a 9 $ 



Considering the estimation accuracy, number of terms, equation order, variable 
selection, etc., several different forms of basis functions were investigated using the 
different resolutions (36, 60, and 360) and all engine operating conditions. Table (5) 
describes each of the basis functions selected from many basis functions that were 
examined. 

Note here that the position function f e for estimating the indicated torque is 
different from the previous one used for the in-cylinder pressure estimation. It is 
effectively a normalized motored torque, which can be calculated from the given engine 
geometry, for each individual cylinder as well as summation of all cylinders. 

Coefficient Training - After selecting one of the prescribed basis functions in 
Table (5), the polynomial coefficients were obtained by taking the same procedures, as 
described in Eqs. (7) and (8). Then, the instantaneous value of the measured basis 
variables or their combinations were used to evaluate each of the polynomial equations 
shown in Table (5) to estimated the desired indicated torque. For instance, choosing 
the basis function 3 would result in the following cross-correlation matrix system. 
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The coefficient set in each basis function was computed once for all at each 
engine operating condition for different number of measurement resolutions. Figures 
(7) and (8) represent the estimated indicated torque in comparison with the measured 
indicated torque using the basis function 3 and 36 samplings per crankshaft rotation at a 
certain engine operating point. 

Referring to Figs. (7) and (8), the indicated torque estimations, either for 
individual cylinder or summation of all cylinders, also provide good agreements with the 
calculated indicated torque traces even based on 36 measurement resolutions. 

Figure (9) represents the estimated indicated torque along with the calculated 
values averaged over each engine cycle, which provides another indication of an 
accurate estimation result using the stochastic approach. The same procedure was 
then applied to 60 resolutions and the other cases of basis functions, and their R.M.S. 
errors are plotted in Fig. (10). Note that the errors indicate the average R.M.S. errors 
over all available engine operating conditions. 

TORQUE ESTIMATION USING FREQUENCY ANALYSIS METHOD 
The goal of this method is to show how crankshaft velocity fluctuations can be 
used to estimate the indicated torque produced by the engine. As explained previously, 
processes involved in generation of the torque are strictly periodic if considered in the 
crankshaft angle domain. The periodicity of the processes suggests the use of Fourier 
Transform as a tool to perform the speed deconvolution through the engine-crankshaft 
dynamics. Again, the approach for the present invention is based on the simultaneous 
measurement of crankshaft speed and indicated pressure in the crank angle domain, 
and on the classical method of frequency identification (experimental transfer function). 
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Based on the SISO model previously described in Fig. (1) and Eq. (4), the spatial 

spectra for the indicated torque and crankshaft speed fluctuations can be constructed 

as shown in Figs. (11) and (12). The first few harmonics of the engine firing frequency 

for these two signals contain enough information in order to represent the actual engine 

behavior, as the firing frequency being defined by the following equation where N is the 

number of cylinder, and S is the stroke. 

' N*2 

A f =— (11) 
f S 

The easiest way to evaluate HQX) at each frequency is to calculate the ratio 
between the DFT (Discrete Fourier Transform) of T e ((jA.) and Q((jX). Instead, a more 
accurate approach takes the measurement noise into account and gives the estimation 
of frequency response of a system using the classical frequency domain estimation 
technique for a SISO system. Using the notation proposed by Bendat and Piersol 
results the following. 

Lower bound for the true frequency response: 

H,=^- (12) 

Upper bound for the true frequency response: 

H 2 =^- (13) 

where Gtt and Gqq are the auto-power spectral densities of indicated torque and 
crankshaft speed while G Tn is the cross-power spectral density between these two 
signals. These quantities are defined as follows: 

Indicated torque auto-power spectral density: 

j M 
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Crankshaft speed auto-power spectral density: 



Speed-torque cross-power spectral density: 



7 M ■ ,x 

Aft/' ' 



To obtain a better estimate of the frequency response the arithmetic average of 
and H 2 has been used such that, 

Arithmetic average of Hi and H 2 : 

H 3 = H,+ 2 H2 (17) 
SI The first few harmonics of the engine firing frequency are sufficient to describe 

the engine behavior. Another reason to use only those components within the entire 
spectra results immediately observing the coherence function between the angular 
velocity fluctuations and indicated torque. Coherence is defined as the following: 



W Coherence function: 

0<y 2 Tn {X)<l (19) 

Because the coherence function gives a measure of how input and output of a 
system are related at a given frequency, it is appropriate to use those frequencies in 
which the coherence is close to one in order to avoid errors due to acquisition noise. 
Figure (13) gives an example of coherence function between indicated torque and 
crankshaft speed fluctuations, and confirms that it is appropriate to use only the first few 
harmonics of engine firing frequency to represent the examined process. Substituting 
values of the crankshaft speed DFT, QQX), and frequency response, H 3 (jZ), in Eq.(4) 
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makes it possible to obtain an estimation of indicated torque. However, this calculation 
does not provide enough information on the average component of the torque. 
Nevertheless, it is possible to extract information on the average torque from its 
fluctuating portion. 

Fourier analysis has shown that the first few harmonics of the engine firing 
frequency can fully describe the fluctuating behavior of the indicated torque as shown in 
Figs. (11) and (12). Experimental results also show that a relationship exists between 
this fluctuating component and the average one. In practice, each variable capable of 
g converting the torque fluctuations as a constant is a candidate to represent this 



HI 



relationship. In this study, the value used for this purpose is an estimate of R.M.S., 
obtained from the following relation, 

j M 

T RMSapprox = ~TT^ T (Ai ) (20) 

Il where M is the number of harmonics taken into the account. Particularly for the 
h average purpose, the first harmonic is considered in the estimation of the average 
torque as shown in the following equation. 

Figure (14) shows the average torque plotted versus the approximated value of 
the R.M.S. Each point in the graph corresponds to a different operating point for the 
engine, with speed varying from 1000 to 5000 RPM. A relationship that is interesting is 
found to be strictly linear at each operating point, and the best-fitted line obtained with 
the least squares method is shown in Eq. (22) below, 

T average = m ' ^RMSappmx + ^ (22) 
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where m = 0.5854 and b = -34.377. This result allows a very important consideration, 
which is an estimate of both fluctuating and average torque components can be 
obtained from crankshaft speed fluctuations only. Also, Fig. (15) shows an example of 
the results obtained from the engine and dynamometer setup at a certain operating 
condition during the experiments. 

REAL-TIME TORQUE ESTIMATION 
The methodology behind the real-time torque estimation is presented with the 
simulation results. Then, the experimental results of the real-time estimation on the 
y, current engine and dynamometer set up are provided as well. The stochastic 

p estimation approach described previously was implemented in real-time. 

I. 

ill Coefficient Estimation - The cross-correlation functions as well as the coefficient 

y 

Mf set in the basis functions were constructed for each specific cases as well as each 

L engine operating condition. In other words, the coefficient set for each basis function is 

w . 

id 

g valid for one specific case and operating condition for which they are evaluated. 

(it 

S| However, in an actual engine operation, these conditions (engine speed and load) are 

s 

continuously changing. To be able to implement the stochastic estimation technique in 
a real-time basis, the indicated torque is estimated accurately over a wide range of the 
engine operating conditions such as speed and load. The pre-computed coefficient set 
of the selected basis function may be stored as a mapping format so that the indicated 
torque may be estimated based on this pre-stored coefficient map at each instance of 
the engine operation. In another approach, each of the basis function coefficients 
themselves is estimated as another function of the engine operating conditions such as 
speed, load, or spark advance. 
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In order to achieve the coefficient estimation technique properly while eliminating 
the need for a coefficient mapping, another set of estimation functions may be 
established that relate each of the coefficients in a basis function to the engine 
operating conditions. Table (6) describes this set of estimation functions, which may be 
specifically used to estimate the basis function coefficients. Note that these estimation 
functions will be referred as "Sub-Basis Functions." In Table (6), 'rpm' represents the 
mean engine speed in RPM, 'Itq' represents the mean engine load, expressed as the 
intake manifold pressure in kPa, and '0 S ' represents the spark advance timing in crank 
angle degree. 



Table 6. Various Sub-Basis Functions 



Function 
Number 


Sub-Basis Function 


1 


a i + b U ' r P™ + b 2,l ' lt <l 


2 


a, = K + K ' r P m + • + K ' r P m ' lt <i 


3 


a t -b 0l +b h -rpm+b 2l -ltq + b 3l - rpm 2 +b 4l -Itq 2 


4 


a, = b 0 j + b l l • rpm + b 2 l • Itq + b 3l * rpm • Itq 
+ b 4l -rpm 2 +b 5t - Itq 2 


5 


a, =b 0l + b u -rpm+b 2} -ltq + b 3i -rpm-ltq + 
b 4A • rpm 2 + b 5 j - Itq 2 + b 6 j ♦ rpm 2 • Itq 2 


6 


a i = b 0t + b u • rpm + b 2 l • Itq + b 2 l • rpm • Itq 
+ b 4l • rpm 2 + b 5 l • Itq 2 + b 6 } • rpm 2 • Itq 
+ b 7 j - rpm • Itq 2 + b 8 i • rpm 2 • Itq 2 


7 


<*, = K + K ' r P m + K * *? + K ' °s 


8 


<*, =boj+b l r r P m + b 2,,'ltq 

+ K t0 s + ** -rpm-6 t +b s , -ltq-0 s 


9 


a , = K + b u - r P m + K • i*q + K • °s 

+ b 4r rpm*0 9 +b 5r Uq.0 a +b €r 9 2 M 
+ b 7r rpm-9 2 s +b 8>r ltq 2 .$ 2 s 



The coefficients b; shown in Table (6) may be determined by minimizing the root 
mean square error between the trained coefficients and the estimated coefficients as 
shown in Eq. (23) below. 
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s = min\^(a tramedj -a estmated<} ) 2 J (23) 

Then, another set of the cross-correlation matrix system, similar to Eq. (10), may 
be constructed to determine the coefficient set bis. As indicated by the seven basis 
functions shown in Table (5) combined with the nine sub-basis functions shown in Table 
(6) for both 36 and 60 resolutions, the coefficient set may actually be expressed as a 
function of the engine mean speed, mean load, and spark advance using any of the 
sub-basis functions described in Table (6). Figure (16) provides an example where the 
coefficients of basis function 3 are estimated using the sub-basis function 2. Note that 
the coefficient shown in this figure is a? in the basis function 3. 

Referring to Fig. (16), the first sub-figure represents effectively the changes in the 
coefficient a? as a function of mean engine speed and load whereas the second sub- 
figure is simply connecting the lines of the first figure in the order of increasing speed 
and load (from left to right in x-axis). Referring to Fig. (16), the trained coefficient a* 
shows a quasi-linear relationship with the engine speed and load, and as a result, the 
sub-basis function (1 st order linear) is able to produce the estimated coefficient with a 
very good accuracy. 

This kind of quasi-linear characteristics of the coefficient with the engine 
operating conditions may be found in those coefficients of linear terms in basis 
functions. In other words, coefficients in the non-linear terms, such as the cross-terms 
in basis functions, typically do not have this type of convenient quasi-linear 
characteristic with respect to the engine operating conditions. To overcome this 
problem, other sub-basis functions with more complex non-linear terms shown in Table 
(6) may be used for the coefficient estimation. 
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INDICATED TORQUE ESTIMATION 
Simulation In Real-Time - In order to simulate the torque estimation in real-time, 
Simulink™ was used to carry out the simulation tasks on the actual engine experimental 
data set described previously. Figures (17) and (18) represent some of the results 
acquired from the simulation of real-time torque estimation. In this example, the 
estimation was carried out based on the choice of basis function 8, sub-basis function 6, 
and 36 resolutions at 2000 RPM and 30 Ibrft. The other cases of the basis and sub- 
basis functions, number of resolutions, and engine operating conditions were also 
I** investigated using the same approach. Figure (19) shows an example of R.M.S. errors 
II resulted from the estimation of indicated torque at each individual cylinder based on 36 
jl resolutions for all basis and sub-basis functions, averaged over all engine operating 
jjj conditions. In the Fig. (19), the bold straight line represents the variation of R.M.S. 
iiq. errors for which the trained (exact) coefficients were used. 

q As it can be observed in Figs. (17-19), even with the estimated coefficient sets 

in 

O the indicated torque estimation for both individual cylinders and summation of all 

III 

cylinders provide accurate results within an acceptable tolerance. Particularly in Fig. 
(19), it may be easily noticed that R.M.S. errors of the real-time torque estimation 
suddenly increase for the basis functions 5 through 7 while they tend to reduce for those 
basis functions when the trained coefficient are used. This result is due to the fact that 
a higher number of basis function consists of more complex 2 nd order non-linear terms 
inside the equation, which eventually makes the coefficients to become highly non- 
linear with respect to the engine operating conditions. As a result of that, the estimated 
coefficients become less accurate, which then leads the higher value of R.M.S. errors 
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for basis functions 5 to 7 as indicated in Fig. (19). For this reason, basis functions 1 
through 4 were implemented in real-time for the further analysis of torque estimation 
during the actual engine operation. 

Estimation During Actual Engine Operation - In order to achieve the real-time 
estimation properly, the dSPACE AUTOBOX system (DS1003) was used for carrying 
out the necessary computational tasks in real-time during the actual engine operation. 
All the results shown are based on 36 resolutions of measurements per crankshaft 
rotation using the basis function 3 and sub-basis function 2 (refer to Tables 5 and 6). 

The estimation of indicated torque for each individual cylinder was first attempted 
applying the method of stochastic estimation. As described previously, coefficients of 
the torque estimation basis function were first estimated before performing the actual 
estimation of indicated torque. Then, applying these coefficients into the basis function 
at each instance of crankshaft position, speed fluctuation, and acceleration, the desired 
indicated torque was estimated. Figure (20) provides an example of the individual 
cylinder indicated torque, estimated in real-time at 1000 RPM of speed and 10 Ibf-ft of 
load torque, and it is compared to the actual value of indicated torque shown in Fig. 
(21), which was acquired previously at the same engine operating condition. 

Torque may be estimated successfully, even in real-time, using this type of 
estimation approach. The estimated torque has a good agreement with the actual value 
overall. This kind of over estimation around the peak value can be compensated by 
using other basis and sub-basis functions. Using the same basis and sub-basis 
functions as for the individual cylinder torque estimation, the summation of indicated 
torque produced by all four cylinders was also estimated directly. Figure (22) shows an 
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example of torque summation, estimated in real-time while the engine was running at 
1500 RPM of speed and 30 lb f -ft of load torque. Then, Fig. (23) provides a comparison 
with the actual indicated torque, which was acquired previously at the same engine 
operating condition. Again, the two figures indicate that sum of indicated torque for all 
cylinders can be accurately estimated as well as individual cylinder torque. Relatively 
simple estimation models, such as basis function 3 and sub-basis function 2, still 
perform a reasonably accurate estimation while keeping the computational 
requirements minimal during the real-time operation. 

Using the present invention, the engine torque generated by each cylinder in an 
IC engine can be successfully estimated based on the crankshaft angular position and 
111 speed measurements. The Stochastic Analysis and Frequency Analysis techniques 

!■'' € 

•"Fr 

p 1 cover a wide range of operating conditions. Moreover, the torque estimation system 
L- and method are independent of the engine inputs (Air, Fuel, and Spark). The procedure 

ui 

p%- allows estimation of not only the cycle-averaged indicated torque but also the indicated 

m . 

Q torque based on the crank-angle resolution with small estimation errors. Furthermore, 

M 

the procedures show the capability of performing torque estimations based on a low 
sampling resolution, thus reducing the computational requirements, which lends itself to 
the real-time, on-board estimation and control. In summary, the approaches may be 
applied for the event-based control in real-time, while eliminating the need for in-cylinder 
pressure transducers. As a result, it is possible to develop practically implementable 
engine diagnostics and control developments providing the individual cylinder 
combustion control, transmission shift control, cylinder deactivation control, which would 
lead to reduced emissions and lower fuel consumptions. 
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